Many-body properties of quasi-one dimensional Boson gas in a narrow CIR 
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We study strong interaction effects in a one-dimensional (ID) Boson gas in a narrow confinement 
induced resonance (CIR) . In contrast to the zero range potential, the ID two-body interaction in the 
narrow CIR can be written as a polynomial of derivative 5-function interaction on many-body level. 
Using the asymptotic Bethe ansatz, we find that the low energy physics of this many-body problem is 
described by the Tomonaga-Luttinger liquid where the Luttinger parameters are essentially modified 
by an effective finite range parameter v. This parameter drastically alters quantum criticality and 
universal thermodynamics of the gas. In particular, it drives the Tonks-Girardeau (TG) gas from 
non-mutual Fermi statistics to mutual statistics or to a more exclusive super TG gas. This novel 
feature is further discussed in terms of the breathing mode which is experimentally measurable. 
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Introduction. Over the past few decades, experiments 
on ultracold bosonic and fermionic atoms confined to one 
dimensional (ID) geometry [ll-Q have provided a bet- 
ter understanding of significant quantum statistical and 
strong interaction effects in quantum many-body sys- 
tems. The experimental measurements to date are seen 
to be in good agreement with results obtained from ex- 
actly solved models 0, @. In particular, Haller et al. @ 
made an experimental breakthrough in 2009 by realizing 
a stable highly excited gas-like phase - called the super 
Tonks-Girardeau (TG) gas - in the strongly attractive 
regime of bosonic Cesium atoms. Such a gas was pre- 
dicted earlier by Astrakharchik et al. Q on a gas of hard 
rods and by Batchelor et al. [8] from the Bethe ansatz 
(BA) integrable Bose gas with an attractive interaction. 

So far most quasi- ID confined systems have been ex- 
perimentally realized by the wide Feshbach resonance 
with combination of confinement induced resonance 
(CIR). At low energies, the quasi-lD atom-atom inter- 
action can be written as a (5-function potential V(x) = 
giuS(x) with an effective interacting strength g\d — 
co/m. Where cq = 2[a±(a±/a s — C)] -1 is determined 
by the 3D scattering length a s . Here a± is the confin- 
ing length and C = 1.0326 is a constant [9j. Across 
a CIR (a s ~ a±), gid can be tuned from +oo to — oo. 
Consequently, the quantum dynamics, correlations and 
thermodynamics of few quasi-lD systems have been ex- 
perimentally elucidated. 

However, the resonance width has an essential effect 
which affects interacting properties on many-body level. 
It has been recently pointed out [l(| that a 3D narrow 
Feshbach resonance may produce strong interaction ef- 
fects due to the resonance structure of its phase shift. In 
the narrow resonance, many-body properties are not only 
determined by the scattering length a s but also by an ef- 
fective range that introduces a strong energy-dependent 



scattering amplitude of two colliding atoms 111, |l2| . In 
this letter, we find that the ID two-body interaction 
across a narrow CIR leads to significant finite range de- 
pendent interaction effects in many-body properties of 
the model, including universal Tomonnaga-Luttinger liq- 
uid (TLL) physics, thermodynamics, quantum criticality, 
super TG phase and collective modes, etc. 



Two-body scattering property and quasi-lD confined 
system. In narrow resonances, the two-body scattering 
physics is well described by an energy dependent scat- 
tering length a^ 1 (E) — a^ 1 (0) — mr^E, where a s (0) is 
zero energy scattering length and E is scattering energy. 
Here ro denotes the effective range and it could be large 
and negative in the narrow CIRs. By solving the two- 
body problem under a tight harmonic confinement in two 
radial directions, one obtains the effective ID scatter- 
ing amplitude [12j /& = — 1/[1 — i2k/c(k)} with c(k) = 
(c^ 1 +4vfc 2 )~ 1 . Here c^ 1 = aj_[a±/a s (u>±_) - C]/2 and 
v = —a^ro/8. The confining length a± — y/l/(mu±) 
is determined by the trapping frequency ui±. The CIR 
corresponds to the limit cq — > oo. The only difference 
between narrow and wide CIRs is that the interacting 
strength c(k), instead of being a constant, depends on 
the relative momentum of two collision atoms. Under 
this circumstance, the energy dependent scattering am- 
plitude is characterized by the parameter v which has 
a dimension of L 3 . For an infinite wide resonance, we 
have ro —> then v = 0, c(fc) = c and the scattering 
amplitude //. becomes identical to the zero range one [§] . 



Remarkably, we prove that for » / the two-body 
scattering amplitude leads to the following Hamiltonian 
with a polynomial of derivative (5-function interaction on 
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Where m is atom mass and N is the total particle num- 
ber. This potential still has the nature of contact inter- 
action. Therefore, in the knowledge of asymptotic BA 
13], for the domain x\ < ■ ■ ■ < xn, the wave function 
of the model ([IT) can be written as as super positions of 
plane waves [l5| , i.e. 




x 2 , • • • , x N ) = y^ A P 



(2) 



where P runs over all N\ permutations, the coefficients 



FIG. 1. The energy density and chemical potential vs zero 
energy interacting strength r^ 1 for the gas-like phase. The 
curves stand for the values of v — 0, 0.01, 0.02, 0.05 (from 
top to bottom). 



we see that the finite range effect does not play a role 
in weak coupling regime. However, for a strong interact- 
ing regime, v drastically changes the ground energy. For 
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wave functions in other domains can be obtained from 
the interchange symmetry of bosons. The eigenvalue of 
the model ([T]) is given by E — X^Li k j wriere the wave 
numbers {kj} satisfy the following transcendental equa- 
tions 
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Here we denoted 9(p) = 2 arctan[p(c ( ^ 1 
integer for odd N and a half odd integer for even N. 
In this work, we only consider the repulsive scattering 
branch corresponding to the solutions with all kj being 
real. 

Ground state and low energy physics. In the thermo- 
dynamic limit N, L — > oo with finite particle density 
n = N/L, we can define a quasi momentum distribution 
Pi,. At the ground state, is non-zero only in a finite 
interval [— Q, Q] where Q is the "Fermi momentum" of 
the ID interacting bosons. From equation ([3]) we have 
p k satisfying 



This means that the effective finite range v can drive the 
the TG gas with non- mutual Fermi statistics into the 
gas-like phase with mutual statistics either less or more 
exclusive than the free Fermion statistics, i.e. super TG 
gas. Fig. Q]shows strong interaction effects in e(7o, v) and 
p, at different values of v. One sees that the energy and 
chemical potential are significantly reduced as v increases 
in the TG and super-TG gas-like phases (in accordance 
with the wide CIR, we here refer the TG and super-TG 
regimes for 70 — > 00 and 70 — > — 00). It indicates that 
the zero range repulsive interaction strength is reduced 
by the positive v at the collision. 

The elementary excitations of the model ([1} are char- 
acterized by the quantum numbers {Ij} for the BA equa- 
tions ([31), including moving a particle close to the right 
or left Fermi points out sides the Fermi sea and moving 
a particle from the left Fermi point to the right [J] . Here 
we find that the elementary spectra are phonon like in 
the long wave length limit and have the same sound ve- 
locity. The sound velocity v c can be obtained from the 
BA equations © by 
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, and the particle 
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density n = Pkdk. Thus the ground state energy 
Eq and chemical potential p are given by 

f Q 

E =L k 2 p k dk = Ln 3 e(7o,w) (5) 
J-Q 

»=^N= n l 3e ~ 70 ^ + 3 ^J (6) 



where f k satisfies the equation f k = [ J_gK(k—q)f q dq+ 
K(Q—k)] I (27r) and uj p refers to the energy of the elmen- 
tary particle-hole excitations. 

On the other hand, the low energy physics of a wide 
class of ID interacting models can be captured by the 
TLL described by the following effective Hamiltonian 

13: 



H eff = 7T / dx[vj{d x cj>{x)f + v N (d x 9(x)) 2 ] (8) 



where 70 = cq/u, v — vn 3 are dimensionless counterpart 
of cq and v (v < cq). By analysis of the BA equation (Q}, 
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where (j>(x) and 9{x) correspond to phase and density 
fluctuations. For the Galilean invariance, vj — v s K — 
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FIG. 2. Sound velocity and Luttinger parameter K vs r 1 
for the values of v = 0, 0.005, 0.01, 0.02 (solid, dashed, 
dash-dotted and dotted lines). In figure (a), all lines stand 
for result of v a — ,Ji>Jvn whereas the open circles present the 
v c from the elementary excitation spectrums J7|. 



nir/m and vn — v s /K are phase and density stiffness 
respectively. Here K is the Luttinger parameter which 
determines the critical exponents of correlation functions. 
The Hamiltonian (J5J also has a linear spectrum uj p = v s p 
with v s — ^Jvjvn- The velocity = can be ob- 
tained from our exact solution ©. In Fig. [5J we ver- 
ify from v s = v c that the effective model © indeed 
captures the collective low energy physics of model ([1}. 
This conclusion will be further identified from the uni- 
versal leading order finite temperature corrections to the 
free energy. The figure b in Fig. [2] shows influences of 
the effective finite range v on the Luttinger parameter 
K = v s jv N . 

Universal thermodynamics. At finite temperatures, 
the true physical states can be determined from the min- 
imization of the Gibbs free energy. Following Yang- Yang 
method [Hj], the equation of state can be obtained from 
the pressure 



rp I- + 00 

P=— ln(l + e - ek/T )dk 
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where the "dressed energy" e k satisfies the so-called ther- 
modynamic BA (TBA) equation 

rp (-+00 

e k = k 2 - f i-— / K(k-q)hx{l + e- e "/ T )dq.{10) 

This equation provides an analytical way to access full 
thermodynamics of the many-body system ([I]). 

At high temperatures, the fugacity z = e M / T <C 1, 
thus we can perform high temperature expansions from 
the pressure (O, namely, 

T 3/2 

P = p(o) + Ab 2 z 2 + 0(z 3 ) (11) 
V2tt 



where S(k) = arctan[2fc/c(fc)] is the two-body scattering 
phase shift. The noninteracting part is given by P^ = 



FIG. 3. Equation of state of the TG gas (co — , oo) is plot- 
ted by numerically solving equation (|10[) . Left panel: pres- 
sure vs temperature at fixed density for v = (solid line), 
v = 0.003 (dash line). The open circles shows the result from 
analytical result (|13|l . Right panel: entropy in the T — fi 
plane. Where V denotes the vacuum and QCR stands quan- 
tum critical regime. The right solid line indicates a crossover 
temperature at which the linear-temperature-dependent en- 
tropy breaks down. 



— y~ f^°° ln(l — ze~ k / T )dk. The second viral coefficient 
A6 2 in (fT2|) fully agrees with the result obtained from the 
pure two-body physics [l^ . In fact, the second term in 
(fTT|) contains subtle corrections from quantum statistical 
and dynamical interactions as well as the finite range 
effect to the pressure of the Boltzmann gas. 

Furthermore, we may accurately calculated the ther- 
modynamics of the model from (1101) in analytic fash- 
ion using the polylog function. For capturing quantum 
criticality of finite range v, we consider the Tonks limit 
Co — > 00. Following the approach [ljij], we obtain an ana- 
lytic form of the equation of state 
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which is valid for finite v. This analytical equation of 
state is very convenient for experimentalists to evalu- 
ate various thermodynamical quantities without involv- 
ing numerical calculation of the integral equation (|10[) . 
From the figure (a) in Fig. [3j we see that a good agree- 
ment between analytical result (|13j) and the numerical 
result from the TBA equation ([T0|) . The significance of 
this result is that the equation of state (fT3")) allows the 
exploration of TLL physics and quantum criticality in a 
highly precise manner. 

At low temperatures, from (|13j) we further find the 
universal field theory of TTL, i.e., up to v 2 , 
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For fixed density n, the free energy can be written as 

ttT 2 

F(n,T)=E -P(n)— + O(T 4 ) (15) 

where E is the ground state energy and j3(n) is given as: 

/9(n) = — (l + 167r 2 fi-— 7r 4 i; a y (16) 

Thus we have obtained the Gaussian conformal field the- 
ory with the sound velocity v c — l//3(n) and central 
charge C = 1. We further prove that this sound veloc- 
ity does coincide with the calculation from the excitation 
spectrum (|7"[). 

Quantum Criticality. By analyzing the TBA equation 
(|10[) in zero temperature limit, we see that the quantum 
phase transition from vacuum to TLL occurs as the chem- 
ical potential is varied through the critical point /i c = 0. 
Near the quantum critical point, this many-body sys- 
tem is expected to show universal scaling behavior in the 
thermodynamic quantities due to the collective nature of 
many-body effect. Taking cq — > oo, the quantum critical 
phenomena of the model reveal subtle dependence of fi- 
nite range parameter v. The universal scaling functions 
of density and compressibility k near this critical point 
can be obtained from the equation of state (|13j) . i.e. 

n(fi, f) - no(fi, f) = fi+ 1 -^(^p) (17) 
«(£, f ) - «o(A, f) = f i+^QfB^p.) (18) 

\ J_ vz / 

where F{x) = -Lii (-e x )/(2v 1 / 3 ^), Q(x) = 

-0Lu_i (-e x ), f = TV/ 3 , ft = fiv 2 / 3 with the back- 
ground no = kq = read off the critical exponents 
z = 2, v = 1. This scaling form looks similar to the 
one of the Lieb-Liniger model 19]. However, in con- 
trast to the scaling behavior of the hardcore Bose gas 
and the Lieb-Liniger model 0, [20| in the wide CIR, 
here the scalings are rescaled by the effective finite range 
parameter v. Fig. |3Jd shows the universal critical behav- 
ior of entropy density for the Cq — > oo limit. The well 
pronounced fan-shape like phase diagram of criticality 
indicates a universal crossover temperature line that sep- 
arates the relativistic TLL from the non-relativistic free 
fcrmion theory. 

Experimental signatures. The oscillation of fundamen- 
tal modes in a shallow harmonic trap (often call breath- 
ing mode) is of particular interest in realistic experi- 
ments. The breathing mode frequency is sensitive to 
various interaction regimes. Therefore, it can be used 
to precisely determine different quantum phases in cur- 
rent experiments. Despite of the fact that for a narrow 
resonance one always has negative effective range such 
that v > 0, it is also possible to create a large positive 
effective range v < like the case of dipolar interaction 
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FIG. 4. Breathing mode frequency across a narrow CIR. Left 
panel: u) as a function of v in the Tonk limit Co — > oo where 
solid line is from numerical calculation while open circles refer 
to small v expansion. Right panel: w as a function of zero en- 
ergy scattering strength in the super-TG regime where solid, 
dash-dotted, dashed lines refer to v — 0, 0.021, — 0.021. 

induced resonance [2l| . We will show that the finite ef- 
fective range v > or v < have dramatic influences on 
the breathing mode frequency oj across a narrow CIR. 

Suppose the ID system is confined in a shallow har- 
monic trap V ex t(x) — ^mu) 2 x 2 along axial direction. 
The breathing mode frequency u) from sum rule is 
given by w 2 /cj^ = —2(x 2 )/(d(x 2 )/duj 2 ), where (x 2 ) — 
J x 2 p(x)dx/N and p(x) is the density distribution along 
the trap. Here p{x) can be obtained from the BA equa- 
tion Q within the local density approximation. Follow- 
ing the convention used in experiment [6], we prefer the 
quantity A 2 = (aid/i x ) 2 to character the zero energy in- 
teraction strength and v = v/£ x to character the effective 
range effect in the trap. In these equations, we have de- 
fined £ x = a x /VN , aid = -2/cq and a x = y/l/{mu x ). 
In Fig |4}d, we plot the ratio of oj 2 /ui 2 as a function of A 2 
in the super-TG regime with |co| 3> 1 at zero tempera- 
ture. One can see that both the value and position of the 
peak change dramatically for the values of v ~ 0, ±0.02. 
This signature is capable of being measured in sufficiently 
narrow resonances. We also note that although for a 
given resonance the effective range ro is fixed, one can 
tune the value of v in a wide range by changing the con- 
fining length a x . Where v oc 1 /a x . 

Moreover, in the "Tonk limit" cq — > oo, by tak- 
ing small v expansions with the equation (|14j) at zero 
temperature, we obtain ui 2 /ui 2 ~ 4/(1 + av + bv 2 ), 
where a = 2048V^/(35tt 2 ) and b = (17827425tt 2 - 
153092096)/(91875tt 4 ). We compare this asymptotic ex- 
pression with numerical result in Fig f4] and see an excel- 
lent agreement between the two results for small values 
of v. We find a peak structure on the negative side of v 
around v ~ — 0.1. This signature indicates that in the 
Co — > oo limit, a negative v gives an effective super-TG 
gas like phase at the ground state. 

Summary. In conclusion, using exact BA solution we 
have studied significantly different strong interaction ef- 
fects in many-body properties of the quasi-lD Boson gas 
across the narrow CIR. The universal TLL physics, equa- 



tion of state, quantum criticality and high temperature 
thermodynamics have been obtained in terms of the ef- 
fective finite range parameter v. It turns out that this 
parameter essentially changes the quantum statistics of 
the TG gas from free Fermi statistics to mutual statis- 
tics or to the more exclusive super TG gas. Measuring 
the breathing mode enables us to capture these striking 
strong interaction effects. 
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